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Abstract 

It is shown that the only functionals, within a natural class, which 
are monotonic in time for all solutions of the vacuum Einstein equations 
admitting a smooth "piece" of conformal null infinity ff, are those depend- 
ing on the metric only through a specific combination of the Bondi 'mass 
aspect' and other next-to-leading order terms in the metric. Under the 
extra condition of passive BMS invariance, the unique such functional (up 
to a multiplicative factor) is the Trautman-Bondi energy. It is also shown 
that this energy remains well-defined for a wide class of 'poly homogeneous' 
metrics. 
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1 Introduction 



Consider a Lagrangian theory of fields 4> A defined on a manifold M with a La- 
grange function density 

£ = £[0 A , ...,«9 Ml ...^ fe A ] , (1.1) 

for some k £ N, where <9 M denotes partial differentiation with respect to x M . Sup- 
pose further that there exists a function t on M such that M can be decomposed 
as R x S, where £ = {t = 0} is a hypersurface in M and the vector d/dt is 
tangent to the R factor. The proof of the Noether theorem, as presented e.g. 
in 



41| , Section 10.1], shows that the vector density 



has vanishing divergence, E x ^ = 0, when the fields (j>A are sufficiently smooth 
and satisfy the variational equations associated with a sufficiently smooth £ (c/. 



also (This is in any case easily seen by calculating the divergence of the 



right-hand-side of eq. (|L~^).) Here (j) A , ai ... ae = d ai . . .d ae 4> A , and X^d^ = d t . 



In first order theories, that is theories in which £ depends only upon (p A and its 
first derivatives, it is customary to define the total energy associated with the 
hypersurface £ by the formula 



E(E) = / E x dS x , (1.3) 



with dS\ = d\jdx° A . . . dx 3 , where j denotes contraction^. By extrapolation 
one can also use ( |1.3| ) to define an "energy" for higher order theories. Because 
of its origin, the right-hand-side of eq. ( |1.3|) will be called the Noether energy of 
S, associated with a Lagrange function £ and with the vector field X. Now it is 
well known that the addition to £ of a functional of the form 

d x (Y x [<f) A , «9 Q A , ...,d ai ... d^^]) , (1.4) 

where k is as in ( |1.1|) , does not affect the field equations^]. We show in Appendix 
[E] that such a change of the Lagrange function will change E(Y0) by a boundary 
integral: 

£(£) — > J E(S)=,E(S)+ / AE^dS^x, (1.5) 

JdT, 

1 03X23X23 

We use the conventions that do j dx A ... dx = dx A dx A dx , d\ j da; A da; A dx = 
dx 2 A da; 3 , etc. 

2 Here we adopt the standard point of view, that the field equations are obtained by requiring 
the action to be stationary with respect to all compactly supported variations (c/. e.g. |33| for 
a discussion of problems that might arise when this requirement is not enforced). 
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where S aj3 = d a j dp j dx° A . . . A dx 3 , with AE^ X given by eq. flETEp . If <9£ is 
a "sphere at infinity" the integral over <9E has of course to be understood by a 
limiting process. Unless the boundary conditions at <9£ force all such boundary 
integrals to give a zero contribution, if one wants to define energy using this 
framework one has to have a criterion for choosing a "best" functional, within 
the class of all functionals obtainable in this way. As discussed in more detail in 
Section ^|, the vanishing of such boundary integrals will not occur in several cases 
of interest. 

Now the concept of energy plays a most important role in the context of fields 
which are asymptotically flat in light-like directions. An appropriate mathemat- 
ical framework here is that of spacelike hypersurfaces which intersect the future 
null infinity Sf + in a compact cross-section K. For such field configurations it is 
widely accepted that the "correct" definition of energy of a gravitating system 
is that given by Freud [TP]], Trautman [H], 4"5|, Bondi et al. [|7|, and Sachs 



which henceforth will be called the Trautman-Bondi (TB) energy. (Because of 
the difficulty of accessing Refs. [53,^3] we have included an appendix (Appendix 
|A|) which describes those results of 



which are related to the problem at 

should make 
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hand. This appendix, together with the date of publication of [45 



it clear why we are convinced that the name of Trautman should be associated 
with the notion of mass in the radiating regime in general relativity.) There have 
been various attempts to exhibit a privileged role of that expression as compared 
with many alternative ones ( Hi~H3|, |6|, [1C|. [22|. |36| . |37] . |39| , |42fl , to quote a few), but 
the papers known to us have failed, for reasons sometimes closely related to the 
ones described above, to give a completely unambiguous prescription about how 
to define energy at 3 ■ (We make some more comments about that in Section ^|, 
cf. also |4"B|| .) In this paper we wish to point out that the TB energy is, up to 
a multiplicative constant aeR, the only functional of the gravitational field, in 
a certain natural class of functionals, which is monotonic in time for all vacuum 
field configurations which admit a (piece of) a smooth null infinity JT 1- . 

We shall also consider a second, somewhat larger, class of functionals, which 
contains Hamiltonians that arise in an appropriate symplectic framework. (It 
will be seen below that the functionals one obtains from the integrals ( |1.3| ) are 
quadratic polynomials of the appropriate Bondi functions and their derivatives; 
there is no reason for the Hamiltonians to satisfy this restriction.) In that larger 
class we describe all monotone functionals and then among these the further 
requirement of passive super-translation invariance also leads to the TB energy 
as the unique expression. The symplectic framework which is appropriate in the 
context of radiating fields will be described elsewhere. 

It is natural to ask why the Newman-Penrose constants of motion 
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or 

the logarithmic constants of motion of |14|, do not occur in the conclusions of 
Theorem |4.1| . These quantities are excluded by the hypothesis that the boundary 
integrand H al3 which appears in the integrals we consider depends on the coordi- 
nates only through the fields. The Newman-Penrose constants could be obtained 
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as integrals of the form ( |1.2| ) (cf. e.g. |47|]) if explicit r 2 factors were allowed in 
H a P. Similarly logarithmic constants could occur as integrals of the form (|1.2|) if 
explicit 1/ In r or r +l In - - 7 r factors were allowed there. 

This paper is organized as follows: In Section § we review some results about 
"energy expressions" in general relativity, and comment on non-uniqueness of 
those. In Section |3| we find all functionals of the fields induced on J by the 
metric which are monotonic in retarded time, in a large class of natural func- 
tionals. In Section [| we analyse those monotonic functionals which are invariant 
under passive BMS super-translations, and prove our claim about uniqueness of 
the Trautman-Bondi mass. In Section [5] we give a super-translation-invariant 
formula for the Trautman-Bondi momentum, for general cuts of ST. In Section 
H we consider the question of convergence of the Freud superpotential to the 
Trautman-Bondi mass for space-times with a polyhomogeneous £[. Remarkably, 
we find that because of some integral cancellations the Freud integral always 
converges to a "generalized Trautman-Bondi" mass, even for metrics which are 
polyhomogeneous of order 1 (cf. Section |5] for definitions). In Section |7| we briefly 
discuss the potential extensions of our results to a Hamiltonian setting. An ap- 
pendix gives a very short review of Trautman's contribution to the notion of 
energy for radiating metrics, while the remaining four appendices contain some 
technical results needed in the body of the paper. 



Non— uniqueness of the Noether energy for grav- 
itating systems 



As an example of applicability of eq. ( |1.5|) , consider a scalar field in the 
Minkowski space-time, with E = {t = 0}. Assume that satisfies the rather 
strong fall-off conditions 

for (t, x) e E we have d ai . . . d a .<p = o(r~ 2 ), < j < k - 1 , (2.1) 



where k is the integer appearing in ( |1 . 1| ) . In this case the boundary integral in 
( |1 . 5|) will vanish for all smooth Y^'s, as considered in eq. (|1.4j ). This shows that 
the eq. ( |1.3| ) leads to a well-defined notion of energy on this space of fields (what- 
ever the Lagrange function £), as long as the volume integral there converges. 
(That will be the case if, e.g., £ has no linear terms in and its derivatives.) 

Consider, next, the same scalar field in Minkowski space-time, with E being 
a hyperboloid, t = \J 1 + x 2 + y 2 + z 2 . Suppose further that C = V M 0V M 0, so 
that the field equations read 

□0 = 0. (2.2) 

In that case the imposition of the boundary condition ( |2.1| ) does not seem to be of 
interest, as such boundary conditions would be incompatible with the asymptotic 
behaviour of those solutions of eq. fl2.2p which are obtained by evolving compactly 
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supported data on {t = 0}. Thus, even for scalar fields in Minkowski space-time, 
a supplementary condition singling out a preferred E x is needed. 

Now for various field theories on the Minkowski background, including the 
scalar field, one can impose some further conditions on E x which render it 
unique p, I5|. The extension of that analysis to the gravitational field carried on 



in |J also leads to a unique E x (namely the one obtained from the so-called "Ein- 
stein energy-momentum pseudo-tensor"), within the class of objects considered. 
While this is certainly an interesting observation, the hypotheses made in that 
last paper are, however, much more restrictive than is desirable. It seems there- 
fore that for gravitating systems another approach is needed. Let us recall how 
the "Noether charge" formalism described in the Introduction works in that case. 
There exist various variational approaches to general relativity, and depending 
upon the point of view adopted one finds the following: 



1. Let C = a/| det g\R/lQir, where the Ricci scalar is considered as a functional 
of the metric field g^, a symmetric connection T^ 7 , and its first derivatives. 
In that case |§ one finds 

E(E) = -i- / V^ X >iy/\feti\dSx 

= Lj V^v^det^V . (2.3) 

This integral is known as the Komar energy, except that ( |2.3| ) is actually 



half of the expression given by Komar [[30 1 



2. Let C = a/I det R^l/X, where the Ricci tensor is considered as a functional 
of a symmetric connection and its first derivatives, and A is a constant. 
The variational equations for such a theory are the Einstein equations with a 
cosmological constant [EB . The Noether energy gives again J2^| the Komar 



integral (|2.3|) . 

3. Let C = a/| det g\R/l&ii, where the Ricci scalar is considered as a functional 
of the metric field g^ v and its first and second derivatives. In that case 
the value of E(E) is given again by the Komar integral ( |2.3| ) (with a 
"wrong" 1/2 multiplicative factor). 

4. Let C = C(g/j,v, g^a) De the Einstein Lagrange function [|TB|, which is ob- 
tained by adding an appropriate divergence to the Hilbert Lagrange func- 
tion a/| det g\R/16ir. In that case one obtains |19|| 



£(£) = / H^dSfu, , (2.4) 

JdT, 

where H^ v is the "Freud superpotential" for the "Einstein energy-momentum 
pseudotensor" , cf. eq. ( J6.1J) below. 
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Yet another approach, leading to a different energy expression, can be found 



in |35j. 



Consider first initial data for, say vacuum, Einstein equations satisfying the 
usual fall-off conditions at spatial infinity; 

g»u - V = Oir- 1 ), d a9ilv = 0{r~ 2 ) . (2.5) 



In that case both the integrals (|2.3| ) and (|2.4| ) converge. When the integral over 
<9£ in ( |2.3|) is evaluated on a "two-sphere at infinity" in Schwarzschild space- 
time one obtains m/2. On the other hand, under the asymptotic conditions ( |2.5| ) 
the integral ( |2.4| ) coincides with the standard Arnowitt-Deser-Misner (ADM) 
expression for energy, and gives m for that same sphere in Schwarzschild space- 
time. 

Under the asymptotic conditions (|2.5|) , a way to obtain a unique expression 
is given by the symplectic formalism. Namely, one can require that -E(S) be a 
Hamiltonian on an appropriately-defined phase space (cf. e.g. P,[Tl],[29],|3||]). This 



requirement, together with the normalization condition that the Hamiltonian 
vanishes on Minkowski space-time, uniquely singles out the Freud-ADM energy 
as the "correct" global energy for general relativistic initial data sets which sat- 
isfy the "spatial infinity asymptotic flatness conditions". Thus the Hamiltonian 
analysis gives a rather satisfactory way of singling out an energy expression at 
spatial infinity. 

Consider, next, hypersurfaces E which extend to Sf and intersect Sf transver- 
sally. There have been attempts to use symplectic methods to define energy in 
this context (see also Hl5|, p4], |2~6| ) . In particular, the analysis of HHI! shows 



that, under appropriate assumptions, the integral of the time- derivative of the 
TB energy over the retarded time gives a Hamiltonian with respect to a proposed 
symplectic structure. This does not allow one to extract the integrand itself from 
the expression for the Hamiltonian in any unambiguous way, for reasons some- 
what analogous to those described in the Introduction. Moreover in those papers 
one has to assume various decay properties of the fields on J7" for large absolute 
values of the retarded time, which have not been established so far. Finally, as the 
symplectic structure considered in |]-^| has a perhaps less universally accepted 
status than the one considered on standard asymptotically flat hypersurfaces, one 
should perhaps also face the question of uniqueness of the symplectic structures 
involved. For all those reasons we conclude that the framework of fl|-|3| fails to 
demonstrate uniqueness of the TB mass. 



3 Monotonic functionals 

From now on, we shall consider metrics g a p defined on appropriately large sub- 
sets of M 4 , but not necessarily globally defined on IR 4 , and satisfying Einstein's 
equations near S- We shall examine a class of functionals which includes all the 
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cases discussed in section |2], and in particular all functionals differing from the 
Hilbert Lagrangian by a divergence. These functionals have the form 

H[u ,g] = lim / H^[g}dS aP , (3.1) 

P-~*°° JS(t=u +p,p) 

dS a/3 = d a j (9/3 j dx° A . . . A dx 3 , 

where 

H aP [g\{x) = H^(g^(x), d a9fll/ (x), . . . , d ai . . . d ak9flv (x)) . (3.2) 

for some I; 6 N, and H a @ is a twice continuously differentiable function of its 
arguments. Here S(t, p) denotes a sphere r = a/ (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 = p, t = 
x° = r. The metrics g a p will be assumed to satisfy the standard fall-off conditions 
corresponding to asymptotic flatness at null infinity. More precisely, consider a 
space-time (M, g) which admits a conformal completion (which in this section 
we consider to be smooth) in the following sense: there exists a manifold with 
boundary (M, g), a diffeomorphic embedding $ : M — > M \ dM, and a smooth 
function Q on M such that §*(Q~ 2 g) = g. We shall also assume that £1\qm = 0? 
that dQ is nowhere vanishing on dM, and that J7 = dM is diffeomorphic to / x S 2 
where I is an interval (possibly but not necessarily equal to K). By a standard 
construction we can introduce Bondi coordinates near J7" (c/. e.g. or |Ii| ) so 
that we have 

ds 2 = dw 2 - 2e 2/3 dw dr + r 2 /i afe (dx a - U a du)(dx b - U b du) , (3.3) 

r 

with x a = (6 1 , 0). We can introduce quasi-Minkowskian coordinates by setting 

u — t — r, x = r sin 9 cos <fi, y = r sin 9 sin 0, 2 = rcos0. (3.4) 

We shall consider only vacuum metrics; recall that this implies the following 
behavior of h ab , (3, U a and V ( 0,|D, 



hot = h b {l + ^ 2 X cd X^+^ + 0(r-"< 



R 1 h ab h cd XacXbd , . 
^ = -32 ^ + ° (r 



-3n 



_ 1 V bX ab ^ 32N a (v) + P a ( X cd Xcd ) + Sxm^ _ , 

- 2^ 2_+ 16^ + lr j ' 1 j 



1/ = r — 2M(f ) H h 0(r J , 

lor 

Here (u) = (w, x a ) and /i a fcdx a dx b = d9 2 +sm 2 9 d0 2 ; X> a is the covariant derivative 
operator defined by h ao . Indices a, b, etc., take values 2 and 3, and are raised 
and lowered with h ab . The tensor field Xab satisfies the condition 

h ab X ab = 0, (3.6) 
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and no other conditions are imposedQon Xabiy) by the vacuum Einstein equations. 
The functions M and N a satisfy the following equations 

® M _ ^lacibd- • , i-r, -r, -ah 
-TT- — ~o 11 XabXcd + T LJ a D h X , 

Ou 8 4 

f)N a 1 
3^— = -V a M + -e ab V b X - K a , (3.7) 
Ou 4 

K a = ^-x\V c x bc +\x cd V d x a c , 
A = h bd e ac V c V b Xda , 

where x '■= duX- Here e a b = d a ^db-sd 2 fj, where d 2 /i = sin 9 d9 A d<fi = ^e a bdx a A dx b 
is the standard volume form on S 2 . If we fix some Uq G I, then the Einstein 
equations do not impose! any restrictions on the function M(uq, 9, <p) and the 
vector field N a (u , 9, <p) on S 2 . 

Equation (|3.5| ) shows that in the coordinate system (|3.4| ) the metric (|3.3|) is 
of the form 

9^ = 9% + ^ + 9 -^ + 0(r-% (3.8) 

with obvious analogous expansions holding for the various derivatives of g^ when 
an appropriate expansion for the derivatives of h a b is assumed. Here g Q = 
diag(— 1, 1, 1, 1). We can now insert a metric of the form (|3.8| ) into a functional 
of the form ( |3.1| ), and as a further restriction we shall require that H has a finite 
numerical value for all fields g^ v of the type described above. Our hypothesis of 
differentiability of H af3 allows us to Taylor expand H a " to order 2 in terms of 
powers of - g° v , d a (g^ v - g° u ) = d a g^, etc., about g^ = g° u . Note that by 
( |3.2|) the 0, . . . , 0] are constants which are either zero or integrate to 

zero in ( |3.1| ) (otherwise the limit in ( |3.1| ) would be infinite), so that 

H[u, g%] = Vuel. 

The 1/r terms in g^ v and its v derivatives will give at most a quadratic contri- 
bution to H, and the 1/r 2 terms at most a linear one, while the remainder terms 
in the Taylor expansion of H a/3 will contribute nothing in the limit r — ► oo. It 
follows that H can be written in the form 

H = [ h[M } MW, . . . , M« N a , iV a « . . . , N« k \ xab, x2 , ■ ■ • x2?, 6, 0]d 2 /x . 
Js 2 

(3.9) 

Here the addition of a superscript (£) to a quantity denotes the £-th ^-derivative 
of that quantity. The square brackets around the arguments of h are meant to 



3 Note, however, that there may be some restrictions arising from some further global hy- 
potheses if those are made. We emphasize that we do not impose any such global hypotheses 
here. 



s 



emphasize the fact that h is not a function but a local functional of the fields 
which is a differentiable function of M, d a M, . . ., d ai . ..d ak M, M {1 \ d a M^\ 
. . ., d ai . . . d ak M^\ . . ., etc., for some finite number of derivatives in directions 
tangent to S 2 . Note that for functionals (|3.1|) the dependence of H on N a , N a ^\ 
. . ., etc., as well as on derivatives of N a , N a ^\ . . ., etc., in angular directions, 
will be linear because N a comes with a factor r -2 , and we shall henceforth only 
consider such functionals. From a symplectic point of view it turns out to be 
natural not to make the hypothesis that h is a quadratic polynomial of the fields 
and their derivatives, as would be the case for functionals ( 3-l|) , and for this 
reason we shall allow arbitrary differentiable functions h in ( |3.9| ), except for the 
hypothesis of linearity in N a with coefficients of linearity independent of the 
remaining fields. 

Assuming that the metric g^ v is vacuum (at least in a neighborhood of ff) we 
can eliminate the w-derivatives of M and N a in favour of M and w- derivatives of 
Xab, using equations (|3.7|) , so that (|3.9|) can be rewritten as 

H = f h[M, N a , Xab , X { al ■ ■ ■ xi?]dV , (3.10) 

where h is still linear in iV a and its derivatives in directions tangent to S* 2 . By an 
abuse of notation here, we still denote the integrand of H by h, although it will 
in general not coincide with the original h of ( |3.9| ). Before we proceed further we 
need the following result based on the work of Friedrich [ 2(], 21 ] as extended by 
Kannar [271: 



Lemma 3.1 Let M (9, 0) be a smooth function on S , Nq(9, 4>) be a smooth 
vector field on S 2 , and x° ab (0, <P), xlb(°, 0); Xab( d i <t>), ■ ■ ■ ^ any sequence of 
smooth symmetric traceless tensors on S 2 . Then there exists e > and a vacuum 
space-time (M, g) with a smooth conformal completion (in the sense described 
above) (M, g) which has a spherical cut u = u Q of £f + D (u — e, tt ) x 'S' 2 suc h 
that the Bondi functions M(u, 9, (f>), N a {u, 9, (f>) and Xab(u, 9, <f>) satisfy 

lim M(u, 9, 0) = M Q {d, 0) , (3.11) 

lim N a (u, 9, 0) = iV a (#, 0) , (3.12) 

U — >U{) 

VieN limxSK 9, 0) = xlbiO,^). (3.13) 

U — >UQ 

Remarks: 1. Actually the functions Xab{u, 9, 0) can be arbitrarily prescribed 
as functions of (u, 9, 0) on an interval (uq — e, uq] for some appropriate e. The 
above weaker claim is, however, sufficient for our purposes. 

2. The limits lim u ^ Uo in the equations above have been introduced to avoid 
talking about space-times with boundary. 

Proof: For x G [0, 1] and Uq — 1 < u < Uq consider metrics of the form 
g^dx^dx" = -Vx 3 e 2p du 2 + 2e w du dx + h ab {dx a - U a du)(dx b - U b du) , (3.14) 
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We wish to show that we can find e > and a metric g^ v of the form ( 3.14j ) 



defined for x G [0, e], x a G S 2 , u G (uq — e, it ] such that x 2 g tliV satisfies the 
Einstein vacuum equations, for which (|3.11|) - (|3.13 ) hold. We shall construct 



the appropriate solution backwards in u on (it — e, uo) x [0, e] x S 2 by solving 
an asymptotic characteristic initial value problem with data given on the null 
hypersurface Af = {u = u } and on a piece of J7"= {i = 0,«G (tt — 1> u o]} x S* 2 - 
This proceeds as follows: For < x < 1 let 9, 0) be any x-dependent 

family of symmetric non-degenerate tensors on S 2 with h® b (0, 9, 0) = h ab (9, 0) 
(the standard metric on S 2 ), with det h® b (x, 9, 0) = det h ab (9, 0), and with 



dh° ab 



dx 



= X° ab ■ (3-15) 

x=0 



set /i Q &(w = uo,x,9,(f>) = h® b (x,9,(fi). Using the Bondi-van der Burg-Sachs pre- 
scription in the coordinate system (r = 1/x, x a ) [[j],|4(| we can find unique smooth 
functions dj(3(u , x 9, 0), diU a (u , x, 9, 0), d l u V(u , x, 9, 0), d l u h ab (u , x, 9, 0), 
such that eqs. (|3.11| )-( |3.13| ) hold and such that for all iV G N the metric g^ u = 
x 2 9ij,u satisfies = 0({u — u ) N ), whatever the fields j3,V,U,h ab as long as 
those fields and their derivatives assume the boundary values obtained above. 
Indeed, the fields appearing at the right hand side of eqs. ( |3.11| )-( |3.13| ) provide 
precisely the data needed for the construction of a solution of the hierarchy of 
equations obtained by ^-differentiating the Bondi-van der Burg-Sachs equations. 
It follows that any geometric quantities, built out of the metric together with an 
arbitrary finite number of its derivatives, calculated at u = uq for any two such 
metrics will coincide at u = uq. 

Because Xab is symmetric and traceless, h ab can be parametrized as 

/ e 27 cosh(25) sinh(2<5) sin# \ , , 

^ sinh(25)sin# e"^ cosh (25) sin 2 9 J ' ^ > 

where we write 7 = c{v)/r + 0(r~ 3 ), 5 = d{v)/r + 0(r~ 3 ). Let cab' be the 
following tetrad field 

eo = e oo' = —d x , 

ei = e xv = e~ 2 \d u + ^Vx 3 d x + U e d e + U%) , (3.17) 

e2 = e-oi' — —i={ e ~ 1 {'~ cosh<5 + i sinh 5)de + e 7 (sinh S — i cosh<5) cosec 98^) , 
v2 

e3 = eic = (e 2 )* , 

where (e2)* denotes the vector whose coordinate components are complex conju- 
gates of those of €2- From ( 3.16|) - (|3.17|) one can calculate the Newman- Penrose 



quantity a (= r i'oo in the notation of p7H , and = r ooi' in the notation of 
to obtain 

a\ x=0 = c + id . (3.18) 
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In p7| the time sense of eo and e\ is unspecified and e2 is only specified up to 



rotations in the e2-e% plane at points in the intersection of N and Sf (in Kannar's 
notation). Since they are then parallelly propagated in Kannar's treatment, these 
rotations are u and r independent. Note that up to these ambiguities the tetrad 
( |3.17| ) coincides with that used in |27), time-reversed, at x — 0, but will in general 



differ from it at other points. This is irrelevant as far as the value of <j\ x= o is 
concerned because a at x = is calculated using only derivatives of the tetrad 
field tangent to the spheres x = 0, u = constant, so that cr\ x=0 calculated for 
the tetrad (|3.17|) will coincide with that calculated in the tetrad used in f27[] (up 



to a constant factor of modulus 1). The essential point is that c and d give the 
requisite data. 

Now let Xab{u, 9, 0) be an arbitrary one-parameter family of symmetric tensor 
fields on S 2 , with u G (u — 1, u ] such that d l u x a b(u Q , 9, <fi) = Xab(@i 0)> from 
Xab{ u , @, 0) we can calculate a\ x=0 . From d l u g^ u \ u=uo , i = 0, 1,2 (which we have 
already calculated previously) we can determine at u = Uq the remaining initial 
data needed for the Friedrich-Kannar asymptotic initial value problem. The 
existence of an e > (depending on the initial data) and a solution of the vacuum 
Einstein equations g^ u defined for (u, x, 9, 0) 6 (u — e, u ] x [0, e] x S 2 assuming 
those initial data now follows from the main theorem of ||27|| . The property that 
the Bondi functions M, N a , h a b and d x h a b parametrizing the metric g^ u assume 
the desired values on Sf + follows from the uniqueness theorems of [^TJ. □ 

We can now pass to the proof of our main result. 

Theorem 3.2 Consider any functional of the form 

H[g] = [ [h(M, T> a M, V ai . . . V ak M, 
Js 2 

Xabi T^cXabi • • • i T-^ci ■ ■ ■ T-^c k Xabi 

(i) -r, (i) T) 71 ^ 

X-ab ' U cX a b > • • • > U c\ ■ ■ ■ LJ c h X a b ' 



X% b \v cX %\...,V Cl ...V CkX W,x a ) 
+ a a N a + a ab V a N h + ... + a ai ... akb V^ . . . V ak N b ]d 2 ^, (3.19) 



where h is a twice continuously differentiable function of all its arguments, with 
some, say smooth, tensor fields a ai ...a k b on S 2 . If H is monotone non-increasing 
in u for all metrics g which satisfy the vacuum Einstein field equations (with 
M , Xab and N interpreted as Bondi functions appearing in g ), then H can be 
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rewritten as 

H = ! *(M - h ac h bd V a V bXcd , x a )d 2 fx, 

with a differentiable local functional^ whose variational derivative 5^/5 f is 
non-negative. 

Proof: Note first that the tensor fields a ab , . . . , a ai ...a k b, can be set to zero by 
integration by parts and a redefinition of a b , a;& — > with an appropriate o^. 
Calculating the w-derivative of ( |3.19| ) we obtain 

dH f , 5h • A •„ 5/i . 5/i ffc+iK ,9 ,„ 

* = JJlM M + ^ + ^ + ' ' ' + )dV ' <3 ' 20) 

Now 5h/5Xab an d all the terms in ( |3.20| ) except for the last one are independent 
of X^ h+1 K u o) ■ If fih/SXab were non-zero for some k > 1 we could, by Lemma fl.l| , 
find a solution of the vacuum Einstein equations with X^ k+1 K u o) so chosen that 
dH/du > 0, which shows that Sh/Sx^b = for all A; > 1. Setting 

h[M, Xab, x a ] = h(M, V a M } ...,V ai ... V ak M, 

Xabt T^cXabi ■ ■ ■ j T^ci • • • T-^c k Xabi 



(1) _ V y (l) _ p T) y(!) 



. • • • . xS* = o, v cX ^ = o,...,v cl . ..v CkX K ; b > = 0, x a ) , 

we obtain from 5h/5x^ — 0, k > 1, 

#(<?,«) = / (MM, Xafe] + d a iV a )d 2 /i, 
Vs 2 



= / (^ + «^ + #-**)dV (3-21) 



Consider, first, equation Q3.21 ) for solutions of the vacuum Einstein equations 



with Xab\u=u = 0. Equations (|3.21|) and ( |3.7|) then yield 



TST = j s y\^iV a M~- A e ah V b X))d 2 ^ 

= [ ^MV a a a ~^-Xe ab V b a a )d 2 fi. (3.22) 
J S 2 6 12 



To proceed further we need to know a little more about A as defined by (|3.7| ). In 
Appendix [B] we show that the image of the operator Xab £acD c 'D b x ab defined 
on traceless symmetric tensors consists precisely of functions of the form Pijj, 

4 That is, \&(/) is a differentiable function of x a , / and a finite number of its derivatives in 
directions tangent to S 2 . 
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where if) is an arbitrary appropriately differentiable function on S 2 and P is the 
projection operator defined as 



Pip = tp - 



3 



(3.23) 



where the $j form an orthonormal basis of the space of spherical harmonics with 
£ = ($o) an d £ = 1 i = 1, 2, 3). Consequently A runs over all smooth 
functions with no zero or first spherical harmonics as Xab runs over all smooth 
symmetric traceless tensors. This, together with Lemma [3.1| (note that M in 
( 3.22) is arbitrary) shows that dH/du in (|3.22j) will have an arbitrary sign unless 



V a a a = 0, e ab V a a b = ^ a&i , 

8=1 

for some constants It follows that 

1 3 



(3.24) 



(3.25) 



(the fact that the above vector field satisfies (|3.24j) can be checked by a direct 
calculation; the fact that there is only one such vector field is shown in Appendix 
[Bj) . Returning to equation (|3.21|) , we obtain from (|3.25|) and 



dH 

du 



s 2 



s 2 



6M [ 8 XabX A UaUhX 



ab 



\a a K a + -^-Xab 

3 OXab 



d 2 fj, 



I 5M 6 4 dxab 



Xab 



1 Sh 



8 5M 



XabX 



ab 



dV 



Define a new functional ^ by 



(3.26) 



*[/, Xab, x a ] := h[M = f + -V a V bX ab , Xab, x a }. 



(|3.26|) can be rewritten as 



du 8 J S 2 5f JS 2 



I u 1 h 6V 

a c V b X ca + -X ca V b a c + — 
6 4 dXab 



Xab 



dV 

(3.27) 
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dH/du will be non-positive for all Xab if and only if S^/Sf is non-negative, and 
the last integral vanishes, which yields 



-a c V b x ca + - X ca V b a c ) . (3.28) 



&Xab V 6 4' 

We wish to show that d a has to be zero. To do this, fix a smooth / and consider 
Gf[xab, x a ] = J ty[f, Xab, x a ]d 2 p as a functional of Xab- Note that if we endow 
the space of the Xab's with a Sobolev space topology Wk^iS 2 ) with some k large 
enough, then Gf will be a twice differentiate function on that space, and by 
flOgl) we have 

G 'M ■= I i 2 ®cV b x ca + 3 X ac V b a c ) u ab d 2 p , (3.29) 
Vl Js 2 

where G'f [u] denotes the derivative of Gf acting on the symmetric traceless tensor 
v. It follows from Schwarz's Lemma that the second derivative G'i of Gf satisfies 
G'f[r,u\ = G'f[u,T], for all smooth symmetric traceless tensor fields r ab and u ab . 
From ( |3.29p we have 

G" f [T, v] := -i- / (2a c V b r ca + 3r ac V b a c ) u ab d 2 fi. (3.30) 

Letting F := | Yli=i w hh some constants at, we have a a = e ab V b F (cf. 
(EQ5D ). We also have X> a X> fe F = (c/. e.^. @, Lemma 5]), so one gets 

T> a a b = e ab F. Using those identities, by integration by parts one obtains 

Gj[r, u] - G" f [v, r] = -\ [ (a c V b r ca + a b V c r ca + 2Fr ac e b c ) v ah d 2 /i . 

O J S 2 

Since v is arbitrary (traceless, symmetric) we obtain 

TS [a c V b r ca + a b V c r ca + 2Fr ac e b c ] = , (3.31) 

for arbitrary r's. Think of the two-dimensional sphere as a submanifold of 1R 3 . 
By a rotation of the coordinate axes we can always achieve F = A cos#, for some 
constant A. Equation (|3.31| ) at a point po lying on the equator, p = {6 = ir/2, <fi ), 
with a = 6, b = 9 reads 

2AZtyr** = . (3.32) 

Consider the smooth traceless symmetric tensor field T ab dx a dx b = p((d6) 2 — 
sin 2 9(d(p) 2 ) + 2ad9d(j), with p and a - smooth functions on S* 2 , supported near 
the equator, and satisfying p(po) = cr(p ) = 0. Eq. ( |3.32j ) implies 
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for all such functions p, so clearly A = 0, and we finally get 

F = . (3.33) 

Define 

*[/, x a ] = h[f + \v a V bX a \ Xab = 0, x a }. 



( |3~30|) and give 

/" m[M-\v a V bX ab }= ( h[M, xab] = H[g] , 
Js 2 4 Js 2 

which is what had to be established. □ 

4 Super— translation invariance 



Theorem |0| does not quite lead to the Trautman-Bondi mass as a preferred 
quantity in the class of functionals considered in that theorem, as it still contains 
an arbitrary function \1/ of M — ^h ac h bd V gVbXcd and a finite number of its angular 
derivatives. Let us show that the further requirement of passive super-translation 
invariance of H can be used to obtain that desired conclusion. Here the qualifi- 
cation "passive" refers to the fact that we use a different Bondi coordinate system 
but we integrate on the same cut of J7 '■ More precisely, consider a functional H 
as in Theorem 13.2. We can calculate the value of if at a cross-section S 2 for a 



metric g, and compare the result with H calculated on the same cross-section 
of for the same metric with a different Bondi parametrization, differing by a 
(finite, or infinitesimal) BMS super-translation. Let S denote a given cut of 
which in some Bondi coordinate system (u, 6, 0) on S is given by the equation 



u = 0, and set 



H(S) = [ *(M - h ac h bd V a V bXcd ){u = 0, 0, 0)d 2 /i . (4.1) 

Consider another Bondi coordinate system (u,0,(f>) = (u — a(9,(p),9,(p), with 
corresponding functions M, Xali e ^ c - As shown in Appendix (see also 
have 



[AM - X ab \\ab] (% 9, <P) = [AM - x a \\ab + A 2 (A 2 + 2)a] (u, 9, 0) . (4.2) 

The overbar in the left hand side of the last equation denotes the quantity AM — 
X ab \\ab calculated in the barred Bondi frame, using the barred Bondi functions M, 
etc. The requirement that H(S), calculated in the unbarred Bondi coordinate 
system, coincides with H(S), as calculated in the barred Bondi coordinate system, 
gives thus the equation 

Va / ^[M-]v a V bX ab ]d 2 fi= [ %M-\v a V b x ab + \(V a V a + 2)V b V b a]d 2 ii. 
Js 2 4 Js 2 4 4 

(4.3) 
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(It should be emphasized that S is not given by the equation u = 0. We are 
not requiring that the value H(S) of H, calculated on the cut S = {u = 0}, 
coincides with that of H(S). That last condition would be the requirement that 
the value of H does not depend on the cut under consideration, which is of course 
absurd in the radiating regime.) Now, elementary considerations using spherical 
harmonics show that \ = (T> a V a + 2)T>bD b a is an arbitrary function such that 



Px = Xi where P is the projection operator introduced in equation (|3.23|) . If 
we replace a by ta in eq. (|4.3|) , differentiate with respect to t, and set t = 0, we 
obtain thus P|j = 0. It follows that there exist constants w 11 , [i = 0, 1, 2, 3 such 

that |j = (w° + w k nk)/4:TT, nk '■= ^ being an orthogonal (but not orthonormal) 

basis in the space SH 1 of the £ = 1 spherical harmonics. The condition that |j 

be nonnegative gives w° + w k n k > for all nk G S 2 . That will hold if and only 
if w° > \w\, where \w\ = \/5 kl WkWi, so that one may think of u> M as of a future 
timelike vector. We have thus obtained ^(f) = (w° + w k nk)f /47r, and finally 

H = ^-Jjw° + w k n k )(M - \ X ab \\ ab ) d 2 » = ^ Jjw° + w k n k )Md 2 fi . (4.4) 

Equation ( (4.4p has the clear interpretation that H is the Trautman-Bondi mass 
as measured with respect to a frame with time-like four-velocity vector (w°,w l ), 
which can be checked from the transformation properties of Bondi coordinate 
systems under (passive) Lorentz transformations. For completeness we analyze 
that question in Appendix 0. 

The results of this section and Theorem |3.2| imply the following: 



(4.5) 



Theorem 4.1 Let H be a functional of the form 

H[g, u] = lim / H aP (g^ Ul g^ a , g^, ai ...a k )dS a p, 

n^no JS(t,p) 



p^oo 

t — p=U 



where the H a/3 are twice differentiable functions of their arguments. Suppose that 
H is monotonic in u for all vacuum metrics g^ v for which H is finite, provided 
that g^ u satisfies 



9nv = Vtiu + — — h — — — 2 V o[t ) 



h].„(u, 8, 6) h?„,(u, 8, d>) „ 
d ai ... d ak {g, v - ^ I ? - JS^J-lll) = ( r - 2 ), 1 < i < k, 



(4.6) 



ik 



for some C functions h®(u, 9, (j)), a — 1, 2. If H is invariant under passive BMS 
super-translations, then the numerical value of H equals (up to a proportionality 
constant) the Trautman-Bondi mass. 
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Proof: If H is monotonic for all such metrics, then it is monotonic for Bondi- 
Sachs type metrics ( |3.3|) for which a quasi-Minkowskian coordinate system (|3.4j ) 
has been introduced. As discussed at the beginning of Section for such metrics 
( |4.5|) can be written as a quadratic polynomial in the relevant fields, linear in 
N a , so that Theorem |3.2| applies. Now the asymptotic behaviour of the functions 
appearing in the metric (|3.3| ) shows that any quadratic terms in M that could 
possibly survive in the limit r — > oo come with no angular derivatives acting 
on M. The defmiteness of the variational derivative of where \& is given by 



Theorem |3.2| , together with Lemma [3.1|, implies then that \I/ is necessarily linear, 
and the result follows from the argument leading to ( |4.4|) . □ 
Note that the trivial monotone functional, namely H = 0, is contained in the 
result above, the relevant constant of proportionality being zero. 



5 General cuts of [f 

So far we have been considering the TB mass of those cuts of J which are given 
by the equation u = 0. Consider now a cut S of J which, in Bondi coordinates, 
is given by the equation 

S = {u = s(9,<j>)} , 

for some, say smooth, function s on S 2 . Theorem |3.2| , together with the discussion 
of the previous section, suggests that it is natural to define 

m TB (S):=-^- [ (AM~ X ab \\a b )(u = s(e i <f ) ),e, ( f ) )smeded ( f ) (5.1) 

167T J S 2 

V k (S):=-±r- [ (4M- X ab \\ab)(u = s(9,4 > ), 9, 4>)n k sin 9 d9d<P , (5.2) 

107T J S 2 

where n k , k = 1,2,3 denotes the functions sin 9 cos 0, sin 9 sin and cos 9, in that 
order. We have: 

• As observed in Section § (cf. eq. ( [4.4|) ), equation ( |5.1D reduces to the 
standard Trautman-Bondi-Sachs definition when s = 0. 

• It also follows from what is said in the previous section that the quantities 
( |5.1| )-( |5T2] ) are invariant under passive BMS super-translations. 

• Equation ( ^4.4|) together with passive super-translation invariance and the 
discussion of Appendix [D] imply that the quantities (p M ) = (rriTB,p k ) trans- 
form as a Lorentz vector under those boosts which map S into itself. 

• The definitions ( |5.1| )-( |5T2| ) allow us to define a flux of energy-momentum 
through a subset of JT + bounded by two cross-sections thereof. More pre- 
cisely, let i — 1, 2 be two cross-sections of JT + which are graphs over the 
cut u = 0: 

Si = {u = Si(9,(f>)} , 
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and let N C J + be such that dN = s 2 {S 2 ) U s^S 2 ). From the definition 
(|D.6| ) and the relation (|3.7|) we have 

m TB (S 2 ) - m TB {Si) = [ (4M - x a6 ||a 6 ) sin#d#d0 = 



1 



dN 
ab 



, Xa b ,uX a \u^6&ud6d<P , (5.3) 

OZ7T 

which can be thought of as a flux of energy through N. A similar formula 
holds for the space-momentum p k defined by (p. 7] ) : 

p\S 2 )-p k {S l ) = -± r [ (AM - X ab \\ab)n k sin 0d0d0 = 

= -^" / Xa^X^sin^d^. (5.4) 

We note that the existence of a flux formula is a rather trivial property, 
since one can always take the u derivative of any integrand to obtain a 
flux. The interest of the above formulae stems from the fact that Xab,u is 
invariant under (passive) super-translations, so that the fluxes ( |5.3j ) 
also share this property. 



Passive super-translation invariance together with the flux formulae ( |5.3j )- 
( |5.4j ) imply that in a stationary space-time the four-momentum p M defined 
by ( p.l| )-( pl2[ ) is S independent. In particular p M vanishes in Minkowski 
space-time, independently of the cut S. 



6 Polyhomogeneous metrics 

Having established the preferred role played by the Trautman-Bondi mass, it is of 
interest to enquire under what weaker asymptotic conditions one can still obtain a 
definition of mass which is finite and monotonic in u. Recall that in |T4| an ad hoc 
definition of mass was given for all Bondi-type metrics with a "polyhomogeneous 
JT', and that mass was shown there to be monotonic. Similarly it was checked 
in H that for a class of asymptotically flat asymptotically vacuum space-timesQ 
the energy expression defined in ]TU[ converges to an appropriately defined Bondi 
mass. From a field theoretic point of view it is natural to define mass in terms of 
an integral, as considered in Theorem |4.1| , using e.g. the Freud potential, where 
the H a/3 of equation ( |3.1| ) is given by the expression (cf. e.g. fl46| ) 

H» v =W u a X a , (6.1) 

W a = 1 9ap (\ det ^I/V^Xa , (6.2) 

1671-^/1 det g pa \ 

5 The class of metric considered in || includes the metrics polyhomogeneous of order 2 
(see pi and below for definitions). 
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with X M = 5q. Inserting the metric ( |3.14 ) into ( |6.2| ), with X a d a = d u and with 
h a f, parametrized as in (|3.16|) , one obtains via a Sheep |J2| calculation 



UT n X a dS 



fJtU 



u=uo,r=ro 



16tt 



4 - 

-r e 



s2 

2/3 



dU b 



2V + 2re 213 cosh(2 7 ) cosh(25) 
(£/V 7 cosh(2<5) + U* sinh(25) sin(0)) 



dr 

QTJ4> 

+ sin(fl)^— (f/V 27 cosh(25) sin(0) + U e sinh(25)) 



sm(9)d6 d<p 



u=uo,r=ro 



dr 

+r 2 V a U a } 
— I h(r-V) 

8U a 

+2r(l - e w cosh(2 7 ) cosh(25)) - r A e- 2f3 h ab —U 
+r 2 V a U a \ s\n{9)de d(j) . 

J u=un,r=rn 



(6.3) 



dr 



(6.4) 



More precisely, this formula is obtained by "covariantizing" (as described in |TT[ ) 
eq. (|6.2j ) with the following flat background metric rj: 

r]^ v dx^dx v = -du 2 - 2dudr + r 2 (d6 2 + sm 2 (6) d(p 2 ) . 

Eq. ( |6.4j ) is exact; no hypotheses about the asymptotic behaviour of the quantities 
involved have been made. Note that the last term in eq. (|6.4| ) integrates out to 
zero. We shall say that a metric is polyhomogeneous of order k if in the Bondi 
coordinates ( |3.3j ) the functions h ab have a polyhomogeneous expansion (see |L4 
for definitions) in which the lnr terms start at a power r~ k : 



h 



ab 



kb + ^ + 

r 



+ 



h k ab n ln n r hl'^hn^ 



+ 



ab 



+ ... 



Consider first metrics which are polyhomogeneous of order 2. We have then 
7 = 0(r _1 ), 5 = 0(r _1 ) and it follows from the Einstein equations as written 
out e.g. in |[J, Appendix Cfl that (3 = 0( r - 2 ), U a = 0(r' 2 ), dU a /dr = 0(r~ 3 ) 



(6.5) 



and r — V — 0(1). Eqs. (gj) and (gj) then give 

H[u , g] = lim -L / (r - V) sin(0) d6d(j) 

r^oo SlT J 52 



which is the standard Bondi integral. Consider, next, metrics which are polyho- 
mogeneous of order 1. In that case one has 7 = 0(r _1 ln 7V r), 5 = 0(r~ 1 \n N r) 

6 There are unfortunately some misprints in Appendix C of flljj ]: 1) The last term in eq. 
(C4), ir 2 cosec6'(W r i3 +4W 3 ), should be replaced by \r cosec 9(rW 13 + 4W 3 ); 2) in the 8th line 
of Eq. (C6) the factor 4 in front of the term 47172/7 should be replaced by 2. 
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for some N. The Einstein equations imply (see the proof of Prop. 2.1 in p"4| ) 
that p = 0(r~ 2 \n 2N r), U a = 0(r~ 2 ]n N r), dU a /dr = 0(r~ 3 \n N r) and r - V = 
0(\n Nv r) for some Ny- Eqs. ( |3.1|) and (|6.4j ) lead again to (|6.5|) . At first sight 
it appears that the integral at the right hand side of ( |6.5|) might diverge for 
some vacuum metrics which are polyhomogeneous of order 1. However, careful 
study of the leading terms in the Einstein equations shows that those terms in V 
which are linear combinations of In* r are divergences, so that their integral over 
a sphere vanishes. Thus the Freud integral always converges to the monotonic 



mass expression considered in ||14|| . Remarkably, the polyhomogeneous case of 
order k > 1 always has a finite energy. 

Let us mention that for metrics which are polyhomogeneous of order k > 2 the 
Freud integral can be given a Hamiltonian interpretation — this will be discussed 
elsewhere. 



7 Closing remarks 

We have shown that every functional of the fields which is monotonic in time in 
a certain class of functionals for all metrics "having a piece of JT' is proportional 
to the Trautman-Bondi mass. The key ingredient of our proof was the Friedrich- 
Kannar construction of space-times "having a piece of JT' . Now in general the 
space-times we have constructed in the proof above will not have any reasonable 



global properties. For example, in Lemma 3.1 the function M could be chosen 



to be negative. In such a case one expects, from the positive TB mass conjec- 



ture, that the space-time constructed in Lemma |3.1| will have no extension with 
complete Cauchy surfaces. Now the property of having such Cauchy surfaces 
is a starting point of any standard Hamiltonian analysis, and for this reason it 
would be rather useful to have an equivalent of Lemma [3.1| in which well behaved 
space-times are constructed. We expect that a result of that kind can be proved, 
under some mild (yet to be determined) restrictions on the function M (such as 
e.g. positivity), and we are planning to investigate this problem in the future. 

Let us finally mention that using similar ideas to those presented here one 
can prove related results for other field theories, such as e.g. Maxwell theory, or 
for scalar fields. More precisely, for a scalar field one has the following: 

Theorem 7.1 The only functional F, in the class of functionals defined in the 
Introduction, of a scalar field (f) on Minkowski space-time, which is monotonic in 
retarded time for all solutions of the massless linear wave equation, and which 
is a Hamiltonian for the dynamics on a hyperboloid E, is the integral H of the 
standard energy-momentum tensor over E. 



To prove this one uses an equivalent of Lemma |3.1| which, for a scalar field 



on Minkowski space-time, can be easily modified to obtain globally defined so- 
lutions. The question of how to define a symplectic structure for dynamics on 



20 



hyperboloids will be discussed elsewhere | I2| . The requirement that the func- 
tional considered is a Hamiltonian leads to the conclusion that F differs from H 
by a boundary integral. Using arguments similar to the ones presented in this pa- 
per (and actually rather simpler, as the corresponding equations on are much 



simpler in the case of a scalar field) one then proves |13| that all the boundary 
integrands, in the case of the scalar field, which have the right monotonicity prop- 
erties, have to integrate out to zero. Minkowski space-time above can be replaced 
by any Lorentzian manifold which has sufficiently regular conformal completions. 

Let us finally mention that one can set up a Hamiltonian framework in which 
some of the problems related to the Ashtekar-Streubel or Ashtekar-Bombelli- 



Reula approaches, listed in section 0, are avoided |L2]|- Unsurprisingly, the 
Hamiltonians one obtains in such a formalism are again not unique, but the 
non-uniqueness can be controlled in a very precise way. The Trautman-Bondi 
mass turns out to be a Hamiltonian, and an appropriate version of the uniqueness 
Theorem |37| proved above can be used to single out the TB mass amongst the 
family of all possible Hamiltonians. 
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A Trautman's definition of mass in the radia- 
tion regime 

In H45|f | Trautman considers gravitational fields for which a coordinate system 
exists in which the metric can be written in the form 

gfiu = Vfiv + O^r- 1 ), g^ p = hp V k p + 0(r~ 2 ) , (A.l) 
(V - h^ p °h pu )k" = 0(r- 2 ) . (A.2) 

Here the functions hp V satisfy = 0(r _1 ), while the null vector field k v is 
defined as follows: Let a be a spacelike hypersurface, and define n M to be a unit 
space-like vector lying in a perpendicular to the sphere r = const, and pointing 
outside it. Trautman defines k u to be n u + t u , where t denotes a unit time-like 
vector normal to a, such that t° > 0. 

Trautman shows that under the conditions ( |A.1| )-( [A~2"| ) the integral at the 
right-hand-side of the equation 

P>] = lw\dS uX (A.3) 



7 The first chapter of J44] is a slightly expanded version of 
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exists^ and is finite because of cancellations among the divergent terms. Here 
ii n is the Freud potential given in Eq. (|6.2|) . Next, Trautman shows that P^c] 
is coordinate independent in the following sense: Let a new coordinate system 
x' v be given by the equations 

x v _> x < v = x v + a " , (A.4) 

with a 1 ' satisfying 

a u = o{r) , a Vtft = b v k^ + 0(r~ 2 ) , (A.5) 

where 

a v = rjufj,^ , b u = 0(r~ l ) , 

and 

a,v tM > = b u ^k p + 0{r- 2 ) , b V)P = 0(r _1 ) . (A.6) 

Those coordinate transformations preserve the boundary conditions introduced 
above. Trautman notices that under those transformations the integrand in 
changes by terms which are 0(r -3 ), so that P M [cr] itself remains unchanged. 



In Section 4 of E5[ Trautman gives the formula for the total energy and 



momentum, which he calls p M , radiated between two hypersurfaces a and a', 

Vil = P^[a} - P^[a'} = J t^dS v , (A.7) 

under the hypothesis that the energy-momentum tensor of matter fields gives no 
contribution on E. Here 

t/ = r^r + 0(r- 3 ) , (A.8) 

where 

Akt = hT(hr - -^rTV) > ( A -9) 

and « is the constant of proportionality between the Einstein tensor and the 
energy-momentum tensor, and it is clear that the integral over E in Eq. ( A.7 ) is 



defined by a limiting process []. He emphasizes that r is invariant with respect to 
the transformations (|A.4 ) and is non-negative by virtue of ( |A.2j ), so that po > 0. 



For our purposes we need to change the definition of fc M given above: we 
require to be a null vector field satisfying 



8 It is clear that S in flA.3| ) is understood as "a boundary of a at infinity" , defined as far as 
integration is concerned by a limiting process. In the section in which he talks about radiating 
fields Trautman does not give a precise definition of what S is. 

9 In the section in which he talks about radiating fields Trautman does not give a precise 
definition of what E is. In a preceding section of p5[ where boundary conditions appropriate 
for spatial infinity are considered he uses the same equation to show that P M is conserved, and 
in that case he defines S as "a time-like "cylindrical" hypersurface at spatial infinity" . 
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1. is normal to the spheres r = {const}, future pointing and outwards- 
pointing; 



2. k^ satisfies the following asymptotic conditions: 

k° - 1 = 0(r- 1 ) , V - — = 0{r- 1 ) . 

r 

(This is compatible with Trautman's definition if one takes a to be the hypersur- 
face {x° = const} in the coordinate system in which ( |A.1[) — ( |A.2|) hold. However, 
the hypersurfaces we consider here are not of this form.) With this modification 
Eq. ( |A.7|) -( |AT8| ) together with positivity of r are the fundamental statement that 
on hypersurfaces which, in modern terminology, "intersect J7" + " the energy can 
only be radiated away. It should be emphasized that this is a more general state- 
ment than that discussed by Bondi et al. and by Sachs four years later [f7|,|40|l, as 
the boundary conditions (|A.1|) -( |A~^) are weaker than those of 0)0- Indeed, con- 



sider a Bondi-Sachs type metric (|3.3| ), with all the functions appearing there sat- 
isfying the fall-off requirements of If quasi-Minkowskian coordinates are 
introduced via the equations (|3.4j) , one finds that Trautman's conditions (|A.1| )- 
( |A.2p hold with k^ = u^. If a is taken to be the null hypersurface {u = uq} (note 



that with our minor modification of the definition of what k^ is, the hypothesis 
that a is spacelike is not needed any more in the above formalism) the four- 
momentum Pfj,[cr} defined by Eq. ( |A.3| ) gives the Bondi mass as defined in 0.|4C 



If a' is taken to be another such null hypersurface, Eq. ( |A.7j ) yields the Bondi 
mass loss formula (integrated in u) . Further, the coordinate transformations ( |A.4j ) 
comprise the BMS "super-translations" : a super-translation given by Eqs. ( |C.1| )- 
]) below corresponds to a transformation ( |A.4| ) with = 0^(0,0) + 0(l/r), 



for some appropriate functions M (6 I ,0), so that in (|A.5| ) vanishes. 

It should be pointed out that, as discussed in Section | above, the fall-off 
conditions ( |A.l| )- (rA~2j ) allow for a large class of metrics with polyhomogeneous 



asymptotics. Last but not least, using the framework of ||5|] reduces the compu- 



tational complexity of the proof of positivity of mass-loss, as compared to several 
other frameworks, e.g. the Bondi-Sachs one. 



B On some operators on S* 2 

Let us denote by A 2 the Laplace-Beltrami operator associated with the standard 
metric on S 2 , A 2 = V a V a . Let SH l denote the space of spherical harmonics of 
degree I (g G SH l ^=^> A 2 g = —1(1 + l)g)- Consider the following sequence 

V° © V° ^ V 1 V 2 V 1 V° © V° . 

Here V° is the space of, say, smooth functions on S 2 , V 1 - that of smooth 
covectors on S 2 , and V 2 - that of symmetric traceless tensors on S 2 . The various 
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mappings above are defined as follows: 



ioi{f,g) = f\\ a + £a b g\\b , 

i 12 (y) = Va{lb + U6 || a - h ab V\\ c , 

*2l(x) = Xa\\b , 
iw(v) = (v a i\ a ,E ab V allb ) , 

where 1 1 is used to denote the covariant derivative with respect to the Levi-Civita 
connection of the standard metric h a b on S 2 . The following equality holds 

i w o i 21 o z 12 o i 01 = (A 2 (A 2 + 2)) © (A 2 (A 2 + 2)) . (B.l) 

Note that we have iio o i 2 i(x) = (x ab \\abi £bc X a b\\ac) Consider the space V° := 
[SH° © SH 1 ^, where 4_ denotes L 2 orthogonality in L 2 (S 2 ) n C°°(S 2 ). Now the 
operator A 2 (A 2 +2) is surjective from V° to V , so that for any \ E V there exists 
/ G such that A 2 (A 2 + 2)/ = A. Consider the tensor field x = «12 ° «oi((/, 0)), 
then shows that x afc ||afe = ^> which establishes surjectivity of the double 

divergence operator, from the space of symmetric traceless tensors to that of 
functions on the sphere which have no zero and first harmonics. Similarly the 
tensor field \ — H2 ioi{(0, <?)) shows that the map V 2 3 Xab —> £C bX ab \\ac ^ 
[SH° © SH 1 ]^ is surjective. 

To justify our claim, that the vector field a given by eq. ( |3.25| ) is the unique 
solution of eq. (|3.24j) , consider the sequence 

V° © V° ^ V 1 V° © V° . 

It is easy to check that 

Ho o i 01 = A 2 © A 2 , 

so if a a \\ b e ab , a a \\ a e (SH ^ then there exist f,gE {SH ^ such that i i(f,g) = 
a, and they are the unique solutions in (SH ) 1 - of the equations: 

A 2 / = a a |, a , A 2 # = a a \\ b e ab . 

Our claim follows immediately from this observation. 

C Super-translations 

As in Appendix || we use the notation f\\ a = V a f, A 2 = V a V a . 

Consider a super-translation which in an appropriate coordinate system on 
SF reduces to a transformation u — ► u — a(8, <ft), for some, say smooth, function 
a on 5" 2 , with the angular coordinate being left invariant. The super-translation 
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can be extended from ^ to a neighborhood thereof in the physical space-time 
using Bondi coordinates (cf. eq. (|3.3f) ). This leads to the following asymptotic 
expansions (see also [[5], p. 119]): 



l alla ~ 2^ ~ 2 « l|a S& + r'W Vi c ) + • • • , (c.i) 



x a = x a + -a l|a I v afe 

r 



u = u — a a" a aii a H - 



^V^-a^ll^'Vl +••• , (C.2) 



1 A 1 

r = r--A 2 a + - 



ab 



,ab , 



ab 



lab , 



No 



X || 6 CK|| a H- — x «|| a& + — x |U a|| a|| 6 - -a" «|| afe - aq| a a n + 

+ ..., (C.3) 



-(A 2 a) 2 -(A 2 «) l|a «|| a 



where f a f )C is the connection defined by the metric h ao . From those formulae we 
obtain the transformation laws for x an d M: 



M{u = u-a(x a ), x a ) 



1 11 

M + 2 X t\\b a W a + lXt a \\ab + lXtu a \\aU\\b 



(u, x a ) 



Xab{U = u - a{x ), X 



Xab - 2a|| ab + h ab A 2 a 



Consider the quantity X|i a &, where | \ab denotes covariant derivatives with respect 
to the transformed coordinates, d a = d a + a^d^. Note that the occurence of u 
derivatives in d a will introduce u derivatives of Xab in the transformation formula 
for this quantity, and one finds that the combination AM — x ab \\ab has a simple 
transformation law with respect to the super-translations: 



[AM - X ab \\ab] (u = u- aid, 0), 0, 0) = [AM - X ab \\ab + A 2 (A 2 + 2) a] (u, 6, 0) . 

The overbar in the left hand side of the last equation denotes the corresponding 
quantity calculated in the new Bondi frame. Note that while the equations ( JU.lj )- 
( |C.3| ) had only an asymptotic character in 1/r, the last three equations are exact; 
in particular no smallness conditions on a have been imposed. 



D Boost— transformations and 

Let A be a boost-transformation with boost parameter v\ by an appropriate choice 
of space-cordinates we can choose it to act along the z axis. In coordinates 
on Minkowski space-time one has 

u 9 „ 9 ^,u. 

« = — u r T. 2}+°~, tan- = e"taii- + - , 

cosh v — sinn v cos U r 2 2 r 
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with remaining unchanged. It follows that on Sf the boost A reduces to the 
transformation 



u 9 9 

u = — — -, tan - = e u tan - , 4> = <p . (D.l) 

cosh v — sinn v cos 9 2 2 

It is natural to interpret ( |D.1| ) as the definition of the action of the Lorentz boost 
A on $ for general space-times admitting a £[. 

Equation ( |L) . 1| ) leads to the following transformation laws 

du = (cosh v — sinh v cos 9)d u , 

dg = u sinh v sin 9d u + (cosh v — sinh v cos 9)dg , 

— sin 9 — d9 _ „. 

sm9 = — — d9= — -, D.2 

cosh v — sinn v cos 9 cosh v — smh v cos 9 

- cosh v cos 9 — sinh 

cosh v — sinh cos 9 

From ( gjg) one obtains the well known statement, that boosts induce conformal 



transformations of "spheres at infinity": if we denote by ip the transformation 
which takes (9, <p) to (9, 0), then 

ip*h ab = ip~ h ab , (D.3) 

with 

ip = cosh v — sinh v cos 9 . 

We note that ip is a linear combination of £ = and £ = 1 spherical harmonics. 
Set 

r = ipr . (D.4) 
The coordinate transformation ( p.l|) , (|D.4j ) preserves the leading order behaviour 



of all the components of the metric ( |3.3| ). It follows from |43|] (compare also fllf ) 
that (|D.1| ), (|D.4|) can be extended to a neighbourhood of $ while preserving the 



Bondi form of the metric ( |3.3| ), the hypersurface u=0 being mapped into the 
hypersurface u — 0. From (p.l|) , (p.4j ) and (|3.3|) at u = one immediately 
obtains 

17 . (D.5) 

so that 

/ Msin^d^d0= / M(coshz/ — sinh v cos 9) sin ^d^d0 . 
is 2 is 2 

It follows that the knowledge of the I = harmonics of M is not sufficient to 
determine the I = harmonics of M. Let us set 

raT£?L=o : = t~ / M| u=o sin0d^d0 , (D.6) 
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p k \u=o:=^[ M\ u=0 n k sin 0d0# , (D.7) 
Air J S 2 

where n k , k = 1,2,3 denotes the functions sin 9 cos <ft, sin 6 sin <p and cos 9, in that 
order. Equations ( |D.1| ), ( p.4| ) and ( |D.5|) also yield 



M cos # sin 8d6 = M(cosh v cos 9 — sinh z/) sin Odd . 
Consequently we obtain the transformation law 

ttitb = m TB cosh v — p z sinh v , (D.8) 

p z = p z cosh v — rriTB sinh v . (D.9) 

As the choice of the axis along which A acts was arbitrary, the set of numbers 
= i m TB,P k ) transforms as a (contravariant) four-vector under the passive 
action of the Lorentz group on JT. It is therefore natural to interprete ttitb as 
the time component, and the p fc, s as space-components of an energy-momentum 
four- vector p^. We use the qualification "passive" above to emphasize the fact 
that such a simple transformation property holds only for those Lorentz trans- 
formations which map a chosen cross-section of 3 into itself. 



E Changes of the Noether charge induced by 
changes of the Lagrange function 

In this Appendix we wish to derive the transformation rule of the "Noether 
charge" ( |1.5|) , when the Lagrange function is changed by the addition of a term 
of the form (JO]), 

C — >C = C + R, R = d x Y x , (E.l) 

with being a smooth function of the fields and their derivatives up to order 
k — 1. Letting Q be an arbitrary domain of R n with smooth boundary and 
compact closure, we have 

/ Y^dS^ = [ Rd n x . 
Jau Jn 

Integration by parts gives 

= I ^a 5 ^ ' (E-2) 
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where 5R/5(p A is the variational derivative of R, for any smooth fields S(f) A . Eq. 
( |E.2|) still holds with Q = R™ if the 5<j) A 's are compactly supported. In that case 
arbitrariness of the S(j) A, s implies 



which expresses the well known fact that the field equations are unchanged by 
the above transformation of the Lagrange function. It follows that 

'an i=0 \o<P ai ... ai j=0 \ a< P» ai ... ai p 1 ...p j J J 

(E.3) 

It is convenient to choose a coordinate system (x^) = (x 1 ^" 1 ) such that dQ is 
given by the equation x 1 = 0, the v a, s, a = 1, . . . , n — 1 being coordinates on dQ. 
Define 

rb A ' m = fh A 

Yai-.-Cie ' 1...1 a\...a t i 

m times 

Integration by parts in (|E.3|) yields 



d n ~ x v = . (E.4) 



fc-l fc-m-l / nyl N 

/ E E (-i) 4 »ax hr^-^+V - 

^ m=0 i=0 \90 ai '... ai 



As the 50 j4 ' m 's are arbitrary we conclude that 

fc-m-l / X 

£ . . . d a A - it" NO. (E.5) 

Let _E A be the Noether current (|1.2|) corresponding to the Lagrange function 
as in ( |E.1| ). For our purposes it is sufficient to consider vector fields X x which 
are transverse to E. We can choose a coordinate system in a neighbourhood of 
S so that £ is given by the equation x l = 0, and moreover X x d\ = d\. From the 
definition of E x and E x we obtain 

fc— 1 fe— m— 1 



^ = ^ + E E R t^ii m S - 9^ 

m=0 i=0 

fc-l fc-m-l , r,y X V 

= ^ 1 + E E (^i-^-)^S 1 -^- 

m=0 »=0 V 0<pax...a,iJ 
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It follows that 



E x dSx - / E x dS 



x 

k— 1 fc— m— 1 



•^m=0 i=0 V 0<Pa x ...aJ 

„ / fc— 1 fc— m— 2 

+ / r -EE ^K 1 * ( E - 6 ) 

"^ 9S \ m=0 i=0 

fc-m-i-2 / \\ 

E c- 1 )^ • • • M ^^'"^ - ^ — xArf ^ • 

£=0 V C,< Pai...a i bi...b e a J J 

The integral over £ in the right hand side of this last equation vanishes by ( |E.5| ), 
which establishes our claim that the Noether charge of S, defined as J s E x dS\, 
changes by a boundary integral under the change ( |E.1|) of the Lagrange function. 
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